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ABSTRACT 

It  is  now  a  well  established  fact  that  galaxies  undergo  significant  morphological  transformation 
during  their  lifetimes,  manifesting  as  an  evolution  along  the  Hubble  sequence  from  the  late  to  the  early 
Hubble  types.  The  physical  processes  commonly  believed  to  be  responsible  for  this  observed  evolution 
trend,  i.e.  the  major  and  minor  mergers,  as  well  as  gas  accretion  under  a  barred  potential,  though 
demonstrated  applicability  to  selected  types  of  galaxies,  on  the  whole  have  failed  to  reproduce  the  most 
important  statistical  and  internal  properties  of  galaxies.  The  secular  evolution  mechanism  reviewed  in 
this  paper  has  the  potential  to  overcome  most  of  the  known  difficulties  of  the  existing  theories  to  provide 
a  natural  and  coherent  explanation  of  the  properties  of  present  day  as  well  as  high-redshift  galaxies. 
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I.  INTRODUCTION 

It  is  a  generally  accepted  view  that  our  understand¬ 
ing  of  the  mechanisms  and  processes  responsible  for  the 
formation  and  evolution  of  galaxies  is  incomplete.  How¬ 
ever,  the  degree  of  this  incompleteness  does  not  come 
into  sharp  focus  until  we  assemble  our  most  compre¬ 
hensive  and  up-to-date  observational  knowledge,  both 
on  the  large-scale  distribution  as  well  as  on  the  inter¬ 
nal  properties  of  galaxies,  and  compare  these  with  the 
predictions  of  our  working  theoretical  models  to  try  to 
make  a  coherent  picture. 

Take  our  own  home  galaxy,  the  Milky  Way,  as  an  ex¬ 
ample.  It  is  a  typical  field  galaxy  of  Sbc  type  (de  Vau- 
couleurs  &  Pence  1978)  in  a  small  group  environment. 
Existing  theories  offer  several  possible  ways  for  the  for¬ 
mation  of  this  type  of  galaxy.  The  scenario  offered  by 
the  earliest  monolithic  collapse  model  (Eggen,  Lynden- 
Bell,  &  Sandage  1962)  is  that  the  Galaxy’s  mass  distri¬ 
bution  acquired  most  of  its  shape  from  the  very  begin¬ 
ning,  i.e.  about  a  Hubble  time  ago,  and  subsequently 
underwent  only  passive  luminosity  evolution  driven  by 
star  formation,  nucleosynthesis  and  element  recycling. 
One  problem  with  this  static  picture  is  that  the  ob¬ 
served  kinematics  of  the  different  age  groups  of  stars 
in  the  Milky  Way  disk  differ  systematically,  manifest¬ 
ing  as  the  well-known  age-velocity  dispersion  relation 
of  the  solar  neighborhood  stars  (Wielen  1977).  In  the 
primordial  collapse  model,  there  is  nothing  which  could 
account  for  this  secular  increase  of  the  velocity  disper¬ 
sion  of  disk  stars. 

Recent  deep  surveys  have  found  that  galaxies  in  the 
general  field  environment  similar  to  that  occupied  by 
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the  Milky  Way  have  undergone  significant  morpholog¬ 
ical  transformation  over  the  cosmic  time,  following  a 
similar  trend  though  not  to  as  dramatic  a  degree  as 
the  cluster  galaxies  we  will  discuss  next.  It  is  found 
that  more  field  galaxies  are  of  earlier  Hubble  types  in 
the  nearby  universe  than  at  the  higher  redshifts  (Lilly 
et  al.  1998).  There  is  also  a  population  of  so-called 
faint-blue  galaxies,  which  are  in  fact  L,  galaxies  having 
luminosities  and  sizes  similar  to  the  Milky  Way,  which 
exist  at  the  intermediate  redshifts  but  which  have  all 
but  disappeared  in  the  nearby  universe  (Ellis  1997). 

Could  this  observed  morphological  evolution  be  due 
to  the  hypothesized  major /minor  merger  events?  The 
thinness  of  the  Milky  Way  disk,  the  lack  of  a  large  pop¬ 
ulation  of  counter  rotating  stars,  as  well  as  the  smooth¬ 
ness  of  the  age-velocity  dispersion  relation  (Figure  1) 
all  argue  against  either  a  major  merger  or  the  accretion 
of  a  satellite  of  significant  mass  over  the  past  Hubble 
time  (Ostriker  1990;  Wyse  2001).  There  is  indeed  a 
known  discontinuity  in  the  Galaxy  age- velocity  disper¬ 
sion  relation  at  about  11  Gyr  ago  (Binney,  Dehnen, 
&  Bertelli  2000:  Gilmore,  Wyse,  &  Norris  2002).  Al¬ 
though  a  merger  has  been  proposed  as  its  origin,  the 
emergence  of  the  spiral  structure  on  the  disk  at  about 
the  same  time  seems  a  more  likely  cause  for  this  dis¬ 
continuity  and  for  the  creation  of  the  thick  disk.  Fur¬ 
thermore,  the  Milky  Way  bulge  stellar  populations  are 
distinctively  different  from  its  known  satellites  such  as 
the  Magellanic  clouds  (Gilmore  2001),  further  argu¬ 
ing  against  the  building  of  a  significant  fraction  of  the 
Bulge  through  satellite  accretion. 

Is  the  bar-driven  gas  accretion  process  (Kormendy 
1982,  1993:  Pfenniger  &  Norman  1990:  Combes  et  al. 
1990;  Pfenniger  &  Friedli  1991;  Friedli  &  Benz  1993; 
Norman.  Sellwood  &  Hasan  1996;  C'ourteau,  de  Jong 
&  Broeils  1996)  responsible  for  Bulge  building?  It  was 
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Fig.  1  . —  The  age-velocity  dispersion  relation  of  the  solar 
neighborhood  stars.  Points:  observational  data  [squares: 
Wielen  (1977);  triangles:  Carlberg  et  al.  (1985)].  Curve: 
fitted  disk  star  heating  due  to  spiral  structure  (Zhang  1999). 


indeed  shown  in  numerical  simulations  that  gas  accretes 
towards  the  center  after  being  shocked  crossing  the 
galactic  bar  potential;  however,  the  abundance  analyses 
of  stars  in  the  thick  disk  and  bulge  of  the  Milky  Way 
indicate  a  high  [ct/Fe]  ratio  (Gilmore  2001),  where  a 
elements  0  and  Mg  are  created  primarily  in  type  II  su¬ 
pernovae  from  massive  young  stars  and  the  Fe  is  gener¬ 
ated  in  type  I  supernovae  from  lower-mass  stars  of  much 
greater  ages.  The  high  [a/Fe]  abundance  ratio  can  arise 
if  the  Bulge  and  thick-disk  stars  all  formed  within  a 
short  time,  thereby  suppressing  enrichment  in  Fe  which 
requires  longer  timescales.  These  timescale  considera¬ 
tions  suggest  that  most  of  the  Bulge  stars  have  formed 
very  early  on  and  probably  not  far  from  their  present  lo¬ 
cations.  This  limits  the  importance  of  gaseous  inflow  in 
building  the  Bulge,  as  continuous  inflows  would  extend 
the  star-forming  epoch  and  so  enable  Fe  enrichment 
from  the  type  I  supernovae  (B.  Waller  2001,  private 
communication;  Jablonka,  Gorgas  &  Goudfrooij  2002). 

We  now  turn  our  attention  to  dense  clusters.  This  is 
in  fact  the  environment  where  the  morphological  trans¬ 
formation  of  galaxies  was  first  indicated  through  the  so- 
called  Butcher-Oemler  (BO)  effect  (Butcher  &  Oemler 
1978a, b).  When  it  was  discovered,  the  BO  effect  re¬ 
ferred  to  a  bluing  of  colors  for  galaxies  in  the  dense 
clusters  at  the  intermediate  redshifts  compared  to  sim¬ 
ilar  density  clusters  in  the  local  universe,  which  contain 
mostly  red  early  type  galaxies.  Recent  HST  observa¬ 
tions  (Couch  et  al.  1994:  Dressier  et  al.  1994)  have 
been  able  to  resolve  the  morphology  of  the  BO  galaxies 
ans  show  that  they  are  mostly  late  type  disks;  therefore 
the  BO  effect  is  now  considered  not  only  a  color  evo¬ 
lution  effect  but  also  a  morphological  transformation 
effect.  Major  mergers  are  not  likely  to  be  responsible 
for  the  observed  morphological  transformation  of  the 


Butcher-Oemler  (BO)  cluster  galaxies,  due  to  the  high¬ 
speed  nature  of  the  encounters  (Dressier  et  al.  1997), 
as  well  as  the  finding  through  numerical  simulations 
that  dissipationless  mergers  between  preexisting  stel¬ 
lar  disks  cannot  account  for  the  kinematics  of  the  early 
type  galaxies,  especially  the  large  ratio  of  the  rotational 
to  random  velocities  (Heyl,  Hernquist,  &  Spergel  1996: 
Cretton  et  al.  2001);  whereas  the  spiral  disks  in  BO 
clusters  are  found  to  have  formed  most  of  their  stars 
before  the  observed  morphological  transformation  had 
taken  place  (Franx  &  van  Dokkum  2001).  Minor  merg¬ 
ers  are  also  unlikely  to  be  the  main  driver  because  there 
does  not  seem  to  be  a  large  reservoir  of  dwarf  spheroidal 
satellites  in  these  clusters  (Trentham  1997)  to  cause  the 
simultaneous  morphological  transformation  of  the  large 
number  of  BO  galaxies. 

The  rapid  morphological  changes  of  galaxies  in  dus¬ 
ters  are  also  not  likely  to  be  produced  by  a  ram  pressure 
gas-stripping  mechanism  alone  (Gunn  &  Gott  1972), 
since  stripping  could  not  lead  to  a  change  of  bulge-to- 
disk  ratio  (Sandage  1983),  whereas  the  morphological 
transformation  of  the  BO  galaxies  from  the  late  type 
disks  to  SOs  and  ellipticals  requires  such  a  change.  The 
so-called  “harassment"  mechanism  had  been  shown  to 
be  effective  in  stripping  away  the  outer  gas  and  trans¬ 
forming  the  small  late-type  disks  into  early  type  dwarf 
galaxies  (Moore  et  al.  1996),  yet  it  was  shown  to  be 
much  less  effective  on  large  disks  (Gnedin  1999). 

Thus  we  see,  for  both  the  field  and  cluster  galaxies, 
the  existing  theories  are  unable  to  provide  satisfactory 
explanations  of  their  formation  and  evolution. 

II.  SECULAR  MORPHOLOGICAL  EVOLU¬ 
TION  OF  GALAXIES 

During  the  past  few  years,  a  new  mechanism  for  the 
secular  evolution  of  galaxies  has  been  proposed  (Zhang 
1996,  1998,  1999)  which  operates  through  large  scale 
coherent  patterns  in  galaxies  such  as  spirals,  bars  or 
other  skewed  mass  distributions.  This  mechanism  has 
been  overlooked  by  past  workers  in  this  field  due  to 
some  subtle  features  of  its  operation. 

(a)  The  Source  and  Sink  of  Outward  Angular 
Momentum  Transport 

Three  decades  have  past  since  the  publication  of 
Lynden-Bell  &  Kalnajs’  seminal  paper  (Lynden-Bell 
&  Kalnajs  1972,  hereafter  LBK)  demonstrating  that 
spiral  density  waves  in  disk  galaxies  can  transport  an¬ 
gular  momentum  (as  well  as  energy)  outward.  Associ¬ 
ated  with  this  outward  angular  momentum  transport  is 
an  expected  secular  redistribution  of  disk  matter,  co¬ 
inciding  with  the  trend  of  the  entropy  evolution  of  a 
self-gravitating  system;  i.e.,  towards  a  more  and  more 
centrally  concentrated  core  together  with  the  build-up 
of  an  extended  outer  envelope.  However,  this  latter  as¬ 
pect  has  rarely  been  discussed,  if  at  all,  in  the  context 
of  the  LBK  theory  since  the  publication  of  their  pa- 
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per.  One  of  the  reasons  for  this  disparity  is  the  equally 
well-known  second  result  from  the  same  paper:  there  is 
no  interaction  between  a  steady  amplitude  spiral  den¬ 
sity  wave  and  the  basic  state  (i.e.  the  axisymmetric 
part)  of  the  galactic  disk  except  at  the  inner  and  outer 
Lindblad  resonances  (ILR  and  OLR).  There  is  thus  an 
apparent  lack  of  consistency  in  the  LBK  theory  since 
it  is  difficult  to  imagine  a  spiral  wave  constantly  trans¬ 
porting  angular  momentum  outward  without  the  ba¬ 
sic  state  mass  distribution  undergoing  a  corresponding 
change.  One  possible  way  out  of  the  apparent  con¬ 
tradiction  is  if  the  spiral  pattern  is  a  transient  phe¬ 
nomenon  as  was  assumed  by  LBK.  The  outward  angu¬ 
lar  momentum  transport  then  leads  to  a  temporary  and 
short-lived  growth  of  a  wave  train  between  the  ILR  and 
OLR:  since  the  wave  has  negative  angular  momentum 
density  inside  corotation  relative  to  the  axisymmetric 
disk,  the  outward  angular  momentum  transport  leads 
to  its  own  spontaneous  growth. 

However,  studies  of  grand  design  galaxies  in  groups 
(Elmegreen  &  Elmegreen  1983,  1989)  suggest  that  most 
spiral  patterns  last  for  at  least  10  revolutions  if  they  are 
triggered  by  interactions.  Recent  N-body  simulations 
have  also  produced  long-lived  spiral  patterns  in  isolated 
stellar  disks  which  lasted  more  than  10  revolutions  with 
essentially  constant  pattern  speed  and  wave  amplitude 
(Figure  2). 

Such  long-lived  spiral  patterns  are  related  to  the 
spontaneously  growing  modes  in  the  galactic  resonant 
cavity  rather  than  wave  trains.  Furthermore,  the  spiral 
modes  cannot  grow  indefinitely.  The  wave  amplitude 
has  to  be  clamped  at  a  finite  value  as  is  observed  in 
real  galaxies.  The  two  conditions,  i.e.,  the  continuous 
outward  angular  momentum  transport  by  a  long-lived 
spiral  pattern  coupled  with  a  finite  wave  amplitude,  im¬ 
ply  that  at  the  quasi-steady  state  of  the  wave  mode  the 
outward-transported  angular  momentum  cannot  come 
from  the  wave  itself  (for  otherwise  the  wave  amplitude 
will  continue  growing  without  bound),  but  has  to  come 
from  the  basic  state  of  the  galactic  disk,  since  the  wave 
and  the  basic  state  are  the  only  two  subsystems  that  we 
divide  the  disk  into.  The  secular  evolution  of  the  mass 
distribution  of  the  galactic  disk  is  thus  an  inevitable 
consequence  of  the  requirement  of  global  angular  mo¬ 
mentum  conservation  and  the  assumption  of  a  quasi¬ 
stationary  spiral  structure;  or,  to  put  it  in  another  way, 
the  globally  self-consistent  quasi-steady  spiral  solution 
is  maintained  at  the  expense  of  a  continuous  secular 
basic  state  evolution.  The  source  and  sink  of  the  angu¬ 
lar  momentum  transported  by  a  quasi-stationary  spiral 
mode  both  reside  in  the  basic  state:  they  are  the  inner 
and  outer  disk,  respectively  (the  dividing  line  between 
the  loading  and  unloading  of  the  angular  momemtum 
is  at  the  corotation  radius,  as  we  will  discuss  below). 
The  energy  and  angular  momentum  exchange  between 
the  wave  and  the  basic  state  thus  serves  as  a  damp¬ 
ing  mechanism  for  the  spontaneously  growing  density 
wave  mode:  since  the  wave  has  negative  energy  and  an¬ 
gular  momentum  density  inside  corotation,  to  receive 


energy  and  angular  momentum  from  the  basic  state  in 
the  inner  disk  limits  the  wave  growth;  similarly,  since 
the  wave  has  positive  energy  and  angular  momentum 
density  outside  corotation,  to  dump  energy  and  angu¬ 
lar  momentum  to  the  basic  state  also  limit  the  wave 
growth  in  the  outer  disk.  In  the  end  the  nonlinearity 
in  this  exchange  process  helps  to  clamp  the  wave  ampli¬ 
tude  at  a  particular  value  which  is  mainly  determined 
by  the  basic  state  properties  (Zhang  1998). 

What  then  is  the  mechanism  through  which  the 
outward-transported  energy  and  angular  momentum 
are  loaded  onto  the  density  wave  in  the  inner  disk, 
as  well  as  unloaded  in  the  outer  disk?  The  relevant 
mechanism  obviously  has  to  involve  the  interaction  of 
the  basic  state  and  the  wave  mode:  Specifically,  it  has 
to  involve  a  dissipative  energy  and  angular  momen¬ 
tum  exchange  between  the  wave  and  the  basic  state; 
i.e.,  the  loading  of  angular  momentum  onto  the  wave 
from  the  basic  state  inside  corotation,  and  unloading 
of  this  angular  momentum  outside  corotation,  with  the 
wave  itself  being  the  carrier  for  the  angular  momen¬ 
tum  transport.  Such  a  mechanism  for  wave/basic  state 
interaction  was  indeed  found  (Zhang  1996,1998,1999), 
and  we  summarize  the  essential  characteristics  of  this 
mechanism  below. 

In  was  first  shown  in  Zhang  (1996)  that  for  a  self- 
sustained  spiral  mode,  the  minimum  of  the  gravita¬ 
tional  potential  of  a  spiral  density  wave  lags  behind 
the  maximum  in  density  in  the  azimuthal  direction  in¬ 
side  corotation,  and  vice  versa  outside  corotation.  The 
phase  shift  between  the  potential  and  density  spirals 
means  that  there  is  a  torque  exerted  by  the  potential 
spiral  on  the  density  spiral,  and,  at  the  quasi-steady 
state  of  the  wave  mode,  a  secular  transfer  of  energy  and 
angular  momentum  between  the  disk  matter  and  the 
density  wave.  The  existence  of  the  phase  shift  between 
the  potential  and  density  spirals  of  a  self-sustained  spi¬ 
ral  mode  is  partly  a  result  of  the  long  range  nature 
of  gravitational  interaction.  It  is  for  this  reason  that 
a  skewed  bar  or  other  skewed  large  scale  patterns  will 
also  possess  a  phase  shift  and  the  associated  collective 
dissipation  through  essentially  the  same  mechanism. 

The  torque  T(r)  applied  by  the  spiral  potential  on 
the  disk  density  in  an  annulus  of  unit  width  can  be 
written  as  (Zhang  1996,  1998) 

p27T 

T(r)  =  d.L/dt  =  r  -E(r  x  VV)zd</> 

Jo 

=  —  7rmrEi (r)Vi (r)  ■  sm(m<j>o(r)),  (1) 

where  E,  V,  Ei,  Vi  are  the  disk  surface  density,  poten¬ 
tial,  the  spiral  perturbation  density  and  spiral  pertur¬ 
bation  potential  in  the  annulus,  respectively;  L  is  the 
angular  momentum  of  the  disk  matter  in  the  annulus, 
(f>Q  is  the  potential-density  phase  shift,  which  is  found 
to  be  positive  inside  corotation  (potential  lags  density) 
and  negative  outside  corotation  (potential  leads  den¬ 
sity),  and  m  is  the  number  of  spiral  arms.  It  can  be 
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Fig.  2. —  Morphological  evolution  of  an  N-body  spiral  mode  in  a  purely  stellar  disk.  Between  the  adjacent  frames  the 
pattern  rotates  about  120°.  From  Zhang  (1998). 
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seen  from  (1)  that  the  torque  T(r)  is  non-zero  only 
when  the  phase  shift  <j>o  is  non-zero.  Furthermore,  the 
contributions  of  both  the  gravitational  and  advective 
torque  couplings  are  included  in  the  single  torque  inte¬ 
gral  given  in  equation  (1)  at  the  quasi-steady  state  of 
the  wave  mode.  The  proof  of  this  fact  is  given  in  the 
Appendices  of  Zhang  (1998,1999). 

The  energy  and  angular  momentum  exchange  be¬ 
tween  the  disk  matter  and  the  density  wave  at  the 
quasi-steady  state  as  indicated  by  equation  (1)  is  achieved 
through  a  temporary  local  gravitational  instability  at 
the  spiral  arms  (Zhang  1996).  The  length  scale  of  this 
instability  at  the  solar  radius  is  calculated  to  be  about 
1  kpc,  comparable  to  the  length  scale  of  the  giant  HI 
and  molecular  complexes  near  the  Galactic  spiral  arm 
region  (Elmegreen  1979). 

In  Figure  3,  we  show  an  image  of  an  N-body  spiral 
mode  from  Zhang  (1996).  There  we  see  the  density  en¬ 
hancement  at  the  inner  edge  of  the  spiral  pattern  inside 
corotation  (the  corotation  radius  rco  =  30  in  this  case), 
reminiscent  of  the  dust  lanes  observed  at  the  leading 
edges  of  the  spiral  arms  of  physical  galaxies  which  sig¬ 
nal  the  presence  of  gaseous  shocks.  However,  here  in  a 
collisionless  particle  disk  of  the  N-body  simulation  we 
have  obtained  the  signature  of  a  shock  wave,  which  is  a 
phenomenon  generally  attributed  to  a  dissipative  sys¬ 
tem.  In  what  follows  we  present  further  evidence  that 
a  spiral  density  wave  is  in  fact  a  propagating  front  of 
collisionless  shock  (Zhang  1996). 
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Fig.  3  . —  Detailed  morphology  of  a  spontaneously-formed 
N-body  spiral  mode,  showing  the  density  maxima  at  the 
leading  edge  of  the  spiral  pattern  for  locations  inside  coro¬ 
tation  (Zhang  1996). 


In  Figure  4,  the  azimuthal  variations  of  the  different 
disk  parameters  are  plotted.  We  see  from  (a)  that  the 
potential  indeed  lags  the  density  for  this  typical  radial 
location  inside  corotation.  Furthermore,  from  (c)  it  is 
seen  that  Toomre’s  Q  parameter  has  a  clear  minimum 
in  the  higher  density  region  of  the  spiral  arms.  This 
sudden  drop  in  Q  signals  the  presence  of  local  gravi¬ 
tational  instability  at  the  arm  region,  (d)  shows  that 
the  velocity  component  perpendicular  to  the  spiral  arm 
suffers  a  sharp  jump  from  supersonic  to  subsonic  (the 
average  sound  velocity  is  about  0.04,  as  shown  in  (b)), 
further  reinforcing  the  impression  of  the  presence  of  a 
shock. 

The  gravitational  instability  and  the  associated  small- 
angle  scattering  of  the  streaming  stars  at  the  arms  of  a 
self-sustained  spiral  wave  is  what  breaks  the  conserva¬ 
tion  of  the  Jacobi  for  a  single  stellar  orbit  in  a  smooth 
and  steady-amplitude  spiral  potential,  or  equivalently 
the  no-wave- basic-state  interaction  conclusion  of  LBK, 
and  this  allows  the  stellar  orbit  to  display  secular  decay 
or  increase. 

Essentially  the  same  result  can  be  obtained  for  a 
star-gas  two-component  disk  (Zhang  1998).  In  the  past, 
discussions  of  secular  evolution  in  galaxies  have  focused 
on  the  accretion  of  gas  under  the  influence  of  a  central 
bar.  This  originates  partly  from  the  mis-conception 
that  “gas  is  dissipative,  whereas  stars  are  not”.  How¬ 
ever,  as  is  well  known,  the  microscopic  viscosity  in  the 
gas  component  is  inadequate  to  support  a  reasonable 
accretion  rate  even  for  proto-stellar  accretion  disks  (see, 
e.g.  Pringle  1981).  Instead,  the  gravitational  viscos¬ 
ity  due  to  the  collective  dissipation  effect  of  the  non- 
axisymmetric  large-scale  structures  has  to  be  responsi¬ 
ble  even  for  the  accretion  of  the  gas  component.  This 
is  because  gravity  does  really  distinguish  whether  the 
underlying  matter  is  made  of  stars  or  gas.  The  past 
star-gas  two-component  N-body  simulations  have  of¬ 
ten  found  that  the  phase  shift  between  the  stellar  and 
gaseous  densities  are  usually  small  (Carlberg  &  Freed¬ 
man  1985),  especially  in  comparison  with  the  phase- 
shifts  of  these  densities  with  respect  to  their  common 
spiral  potential  (Figure  5).  These  phaseshifts  causes 
stars  and  gas  to  both  drift  towards  the  center  as  well 
as  being  heated.  We  will  calculate  these  spiral  induced 
evolution  rates  quantitatively  in  the  next  subsection. 

(b)  Astrophysical  Consequences 

As  a  result  of  the  wave-basic  state  angular  momen¬ 
tum  exchange,  an  average  orbiting  star  in  the  basic 
state  inside  corotation  loses  energy  and  angular  mo¬ 
mentum  to  the  wave  secularly  and  tends  to  spiral  in¬ 
ward  (Figure  6).  Similarly,  a  star  outside  corotation 
gains  energy  and  angular  momentum  from  the  wave 
and  drifts  outward  secularly  (Figure  7).  The  mean  or¬ 
bital  radius  evolution  leads  to  a  corresponding  disk  sur¬ 
face  density  evolution  which  we  plot  in  Figure  8,  where 
the  dashed  line  in  each  frame  indicates  the  surface  den¬ 
sity  at  the  earlier  time  step,  and  the  solid  line  the  later 
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Fig.  4. —  Spiral  gravitational  shock.  Different  frames  show  the  azimuthal  distributions  of  the  following  parameters:  (a) 
Surface  density  (solid  line)  and  negative  potential  (dashed  line),  (b)  Radial  velocity  dispersion,  (c)  Toomre’s  Q  parameter, 
(d)  Velocity  component  perpendicular  to  the  spiral  arm.  The  above  quantities  are  computed  at  a  radius  of  14.5  (From  Zhang 
1996). 


radius 


Fig.  5  . —  Potential  and  density  phase  shift  for  the  stellar  component  (solid)  and  the  gas  component  (dash),  respectively, 
for  the  spiral  mode  in  the  two-component  N-body  simulations  of  Zhang  (1998)  at  time  step  1600. 
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time  step.  Since  the  star  inside  corotation  drifts  in¬ 
ward,  and  outside  corotation  drifts  outward,  at  corota¬ 
tion  the  surface  density  decreases  with  time.  Note  that 
as  a  galaxy  evolves,  the  spiral  pattern  speed  tends  to 
decrease  and  thus  the  corotation  radius  tends  to  move 
outward  (Toomre  1981).  We  also  see  the  trend  of  in¬ 
creasing  disk  central  density  together  with  the  build-up 
of  the  extended  outer  envelope,  as  predicted. 


Fig.  6  . —  Evolution  of  orbit  trajectory  for  a  typical  star- 
inside  corotation.  From  Zhang  (1996). 


Fig.  8  . —  Evolution  of  disk  surface  density.  From  Zhang 
(1996). 


spiral,  and  m  and  <f> o  are  again  the  number  of  spiral 
arms  and  the  potential-density  phase  shift,  respectively. 

To  calculate  the  secular  evolution  rate  for  our  own 
Galaxy,  we  assume  a  two-armed  spiral  pattern  of  20% 
amplitude  and  20°  pitch  angle  (Drimmel  2000),  which 
is  appropriate  for  the  average  Hubble  type  that  our 
Galaxy  had  during  the  past  1010  years  of  evolution. 
This  set  of  values  gives  an  orbital  decay  rate  of  2  kpc 
per  Hubble  time  using  (2);  the  same  set  of  spiral  param¬ 
eters  also  nicely  fitted  the  observed  age- velocity  disper¬ 
sion  curve  (Figure  1)  using  the  velocity  diffusion  equa¬ 
tion  we  will  derive  later,  increasing  the  credibility  of 
this  particular  set  of  choice  of  spiral  parameters. 

Therefore,  a  star  in  the  Sun’s  orbit  will  not  make  it 
all  the  way  in  to  the  inner  Galaxy  in  a  Hubble  time. 
However,  the  corresponding  mass  accretion  rate  across 
any  Galactic  radius  inside  corotation  is  given  by 


Fig.  7  . —  Evolution  of  orbit  trajectory  for  a  typical  star- 
outside  corotation.  From  Zhang  (1996) 


A  quantitative  estimate  of  the  rate  of  orbital  change 
can  be  obtained  from  equation  (1),  which  turns  out  to 
be  (Zhang  1998) 

dv  1 

—  =  ——F2Vq  tan(i)  sin(m^>o),  (2) 

where  F  is  the  fractional  wave  amplitude,  Vg  is  the 
circular  velocity  of  the  star,  i  is  the  pitch  angle  of  the 


dM  dr  _ 

— —  =  27T  r— E. 
dt  dt 


(3) 


where  E  is  the  disk  surface  density.  Using  (3),  and 
a  solar  neighborhood  average  disk  surface  density  of 
60 MQpc~'2  (Bahcall  1984;  Kuijken  &  Gilmore  1989), 
the  mass  accretion  rate  for  the  Galaxy  disk  is  found  to 
be  about  6  x  109Me  per  1010  yr.  A  substantial  frac¬ 
tion  of  the  Galactic  bulge  can  thus  be  built  up  in  a 
Hubble  time.  The  vertical  drift  (or  velocity  dispersion 
increase)  needed  for  a  star  to  truly  become  a  bulge  star 
is  produced  by  the  isotropic  heating  effect  accompany¬ 
ing  accretion,  which  we  will  discuss  a  little  later. 

Observationally,  bulges  are  found  to  be  old,  and  en¬ 
compass  a  wide  range  of  metallicities  (see,  e.g.  Goud- 
frooij,  Gorgas  &  Jablonka  1999),  with  a  clear  radial  gra¬ 
dient  in  both  the  age  and  metallicity  distribution  which 
appears  to  be  the  continuation  of  a  similar  gradient  in 
the  disk  (Courteau  et  al.  1996).  The  secular  evolu¬ 
tion  process  discussed  in  this  work  predicts  that  the 
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abundance  gradient  will  be  enhanced  over  long  period 
of  time  across  the  entire  galaxy,  since  it  is  effectively  an 
inner  disk  contraction  process  which  amounts  to  about 
2  kpc  of  disk  scalelength  reduction  in  a  Hubble  time. 
This,  coupled  with  the  shorter  dynamical  time  scale  in 
the  inner  disk  and  the  resulting  enhanced  element  recy¬ 
cling  rate,  leads  to  an  increase  of  the  abundance  gradi¬ 
ent  with  time.  This  is  in  contrast  with  the  gaseous-bar 
secular  evolution  scenario  (e.g.  Friedli,  Benz  &  Kenni- 
cutt  1994  and  the  references  therein)  which  generally 
predicts  a  flattening  of  the  abundance  gradient  due  to 
gas  inflow  ,  though  certain  amount  of  flattened  gradient 
is  indeed  observed  in  low  galactic  latitudes  of  the  late 
type  barred  galaxies.  Gas  accretion  also  tends  to  create 
a  bluing  of  the  nuclear  region  even  though  early  type 
bulges  are  often  observed  to  be  red.  Furthermore,  while 
gas  accretion  could  be  important  for  the  morphological 
transformation  among  later  Hubble  types,  it  is  insuffi¬ 
cient  to  explain  the  transformation  of  Sb  to  Sa,  and  Sa 
to  SO,  etc..  Bulge  formation  also  cannot  be  solely  due 
to  the  dissolution  of  a  pre-existing  stellar  bar  as  some 
theories  suggest  since  observationally  the  bulge  light  is 
found  to  be  added  onto  the  disk  light,  instead  of  sub- 
stracted  from  it  (Wyse,  Gilmore, &  Franx  1997).  There 
also  appears  to  be  a  deficiency  of  bars  for  galaxies  be¬ 
yond  z=0.5  (Abraham  et  al.  1999)  and  therefore  cannot 
account  for  the  existence  of  some  high-z  spheroidals.  In 
the  more  complete  version  of  our  secular  evolution  sce¬ 
nario,  star  and  gas  co-accrets  to  the  center  of  a  galaxy, 
therefore  the  inner  disk  effectively  contracts  with  time. 
This  naturally  accounts  for  both  the  late  emergence  of 
bars  and  bulges,  as  well  as  the  fact  that  the  material 
making  up  these  bulges  and  bars  appear  to  have  formed 
stars  relatively  early  and  in  rapid  episodes. 

Another  important  consequence  of  spiral-induced 
wave-basic  state  interaction  is  the  secular  heating  of 
the  disk  stars,  believed  to  be  the  main  process  respon¬ 
sible  for  producing  the  age- velocity  dispersion  relation 
of  the  solar  neighborhood  stars  (Figure  1).  As  we  have 
mentioned  above,  the  secular  heating  process  allows  the 
stars  to  gradually  drift  out  of  the  galactic  plane  as  they 
spiral  inward,  and  eventually  become  bulge  stars. 

The  secular  heating  of  the  disk  stars  works  as  fol¬ 
lows.  Since  a  spiral  density  wave  can  only  gain  energy 
and  angular  momentum  in  proportion  to  the  pat¬ 
tern  speed  of  the  wave,  and  a  disk  star  which  moves  on 
a  nearly  circular  orbit  loses  its  orbital  energy  and  an¬ 
gular  momentum  in  proportion  to  Q,  the  circular  speed 
of  the  star,  an  average  star  cannot  lose  the  orbital  en¬ 
ergy  entirely  to  the  wave:  thus,  the  excess  energy  serves 
to  heat  the  star  when  it  crosses  the  spiral  arm.  For 
our  Galaxy,  the  diffusion  coefficient  due  to  the  spiral- 
induced  secular  heating  is  estimated  to  be 

D  =  (H  —  flp)F2v^  tan(i)  sin(m^&)  ~  6.0 (kms~1)2yr~1 , 

(4) 

if  using  the  same  set  of  spiral  parameters  as  used  above 
for  estimating  Bulge  building  (i.e.  a  20°  pitch  angle 
and  20%  amplitude  two  armed  spiral).  This  value  of 


D  fits  very  well  the  age-velocity  dispersion  relation  for 
the  solar  neighborhood  stars  as  can  be  seen  in  Figure 
1.  The  above  expression  for  D  can  be  shown  to  be 
approximately  independent  of  galactic  radius  (Zhang 
1999),  which  would  reproduce  the  observed  isothermal 
distribution  of  the  stellar  and  gaseous  mass  across  the 
Galaxy  (Gilmore,  King  &  van  der  Kruit  1990).  Since 
the  spiral  gravitational  instability  which  mediates  the 
wave-star  energy  and  angular  momentum  exchange  is 
a  local  instability,  the  heating  of  the  disk  stars  is  ap¬ 
proximately  isotropic,  and  all  three  dimensions  of  the 
velocity  dispersion  increases  at  approximately  the  same 
rate  as  is  observed  (Wielen  1977). 

The  velocity  dispersion  of  the  gas  in  the  high  red- 
shift  Damped  La  systems  (DLAs),  which  are  believed 
to  be  the  candidate  primordial  disk  galaxies,  is  found 
to  be  around  10  km/s  (Wolfe  2001).  This  can  be  grad¬ 
ually  increased  to  the  stellar  velocity  dispersion  of  40 
km/s  of  the  thick  disk  stars  of  the  present-day  Milky 
Way  like  galaxy  through  the  above  spiral  heating  mech¬ 
anism.  The  metallicity  of  the  DLA  systems  are  found  to 
have  evolved  little  between  z=2-4  (Prochaska,  Wolfe  & 
Gawiser  2000),  consistent  with  the  fact  that  the  disks 
during  this  period  have  not  formed  prominent  spiral 
patterns  and  thus  the  metal  enrichment  evolution  is 
not  prominent.  Emergence  of  the  spiral  heats  the  disk 
immediately,  as  confirmed  in  the  N-body  simulations. 
This  can  be  a  natural  explanation  for  the  discontinu¬ 
ity  found  in  the  age-velocity  dispersion  relation  of  the 
solar  neighborhood  stars  11  Gyr  ago. 

A  similar  energy  injection  into  the  interstellar  medi¬ 
um  can  serve  as  the  top-level  energy  source  to  power 
the  subsequent  supersonic  turbulence  cascade  (Zhang 
2002;  Zhang  et  al.  2001),  which  naturally  explains  the 
size-linewidth  relation  of  the  interstellar  clouds. 

We  thus  see  that  the  spiral  and  bar-induced  radial 
mass  accretion  process  leads  to  the  building  up  of  the 
bulge,  and  causes  the  Hubble  type  of  a  galaxy  to  evolve 
from  late  to  early.  Such  morphological  transformation 
is  observationally  most  pronounced  in  dense  BO  clus¬ 
ters,  though  its  less  pronounced  counterpart  in  the  field 
has  also  been  observed  (Lilly  et  al.  1998).  The  en¬ 
hanced  mass  accretion  rate  for  cluster  galaxies  is  in 
part  due  to  the  large  amplitude  and  open  spiral  pat¬ 
terns  induced  through  tidal  interactions  of  neighboring 
galaxies,  since  the  effective  evolution  rate  is  propor¬ 
tional  to  wave  amplitude  squared  and  the  spiral  pitch 
angle  squared  (equation  2)  (note  that  the  phase  shift 
(fro  itself  is  approximately  proportional  to  spiral  pitch 
angle).  Preliminary  evidence  has  been  found  that  the 
brightness  (or  amplitude)  of  the  spiral  pattern  is  signif¬ 
icantly  higher  for  cluster  galaxies  compared  to  the  field 
galaxies  of  the  same  rotation  curve  and  size  (Aragon- 
Salamanca  et  al  2002). 

Recently,  a  class  of  so-called  “passive  spiral  galaxies” 
have  been  discovered  (Goto  et  al.  2003  and  the  ref¬ 
erences  therein),  which  have  spiral  morphologies,  but 
do  not  show  any  star  formation  activity.  The  char- 
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acteristic  environment  in  which  these  galaxies  reside 
coincides  with  the  environment  where  galaxy  star  for¬ 
mation  rates  suddenly  decline  and  the  morphological- 
density  relation  turns.  The  existence  of  passive  spiral 
galaxies  suggests  that  a  physical  mechanism  such  as 
secular  evolution  that  works  calmly  is  preferred  to  dy¬ 
namical  origins  such  as  major  mergers  since  these  can 
destroy  spiral  arm  structures.  Passive  spiral  galaxies 
are  likely  to  be  a  galaxy  population  in  transition  be¬ 
tween  red,  elliptical/SO  galaxies  in  low  redshift  clusters 
and  blue,  spiral  galaxies  more  numerous  in  high  red- 
shift  clusters.  They  are  also  likely  to  be  the  slightly 
evolved  conterpart  of  the  so-called  “E+A”  galaxies  in 
clusters  (Belloni  1997  and  the  references  therein)  which 
show  galaxy-wide  enhanced  star-formation  episode  1- 
2  Gyr  prior  to  the  observation  epoch.  The  fact  that 
the  number  density  of  these  passive  galaxies  seems  to 
be  highest  at  about  1-2  virial  radii  from  the  cluster 
center,  where  the  “SFR-Density”  and  “Morphology- 
Radius”  relations  turn  (Nichol  2003:  Goto  2003)  also 
finds  natural  explanation  in  the  secular  evolution  sce¬ 
nario:  for  cluster  radii  much  outward  of  this  critical 
radius,  the  density  of  galaxies  are  insufficient  to  cause 
significant  tidally-induced  spiral  activity  enhancement 
through  neighboring  galaxy  interactions.  Whereas  for 
cluster  radii  much  smaller  than  the  critical  radius,  the 
tidally  enhanced  star-formation  happened  in  the  long 
past,  and  the  subsequence  secular  evolution  has  trans¬ 
formed  the  morphology  of  these  previously  disky  galax¬ 
ies  into  SOs  and  ellipticals. 

Tidal  interactions  among  neighboring  galaxies  have 
also  been  found  to  produce  enhanced  disk  mass  accre¬ 
tion  and  nuclear  activity  for  galaxies  in  small  groups 
(Byrd  et  al.  1986;  Zhang,  Wright,  &  Alexander  1993). 
The  necessity  of  invoking  environmental  effects  to  en¬ 
hance  the  strength  of  the  spiral  structure  does  not  make 
the  process  we  discuss  less  interesting  or  relevant.  Just 
as  in  the  case  of  the  development  of  a  human  being, 
even  though  the  environmental  input  and  nourishment 
are  important,  without  an  innate  mechanism  for  de¬ 
velopment  and  maturation,  the  human  growth  process 
would  have  been  impossible. 

(c)  Formation  of  Coherent  Patterns  in  Disk 
Galaxies  as  an  Example  of  Non-Equilibrium 
Phase  Transitions 

It  is  well  known  that  for  an  isolated  system,  the 
direction  of  entropy  evolution  is  towards  an  increas¬ 
ing  degree  of  macroscopic  uniformity,  corresponding 
to  increasing  entropy.  For  open  systems  at  far-from- 
equilibrium  conditions,  however,  it  often  happens  that 
the  usual  near-equilibrium  thermodynamic  branch  of 
the  solution  becomes  unstable,  and  new  types  of  highly 
organized  spatial-temporal  structures  emerge  sponta¬ 
neously.  Due  to  its  similarity  to  equilibrium  phase  tran¬ 
sitions,  this  kind  of  spontaneous  structure  formation  in 
nonequilibrium  systems  has  been  termed  “nonequilib¬ 
rium  phase  transitions” ,  and  the  structures  thus  formed 


“dissipative  structures”  (Glansdorff  &  Prigogine  1971: 
Nicolis  &  Prigogine  1977)  to  emphasize  the  constructive 
role  of  dissipation  in  the  maintenance  of  these  nonequi¬ 
librium  structures. 

The  large-scale  coherent  patterns  formed  in  open 
and  nonequilibrium  systems  are  functional  as  well  as 
architectural.  One  of  the  important  functions  of  these 
“dissipative  structures”  is  to  greatly  accelerate  the 
speed  of  entropy  evolution  of  these  systems  towards 
reaching  thermodynamic  equilibrium,  or  at  least  re¬ 
ducing  the  degree  of  nonequilibrium.  The  local  highly- 
ordered  structure  (which  has  low  entropy)  maintains  its 
constant  entropy  in  the  meta-stable  state  by  continu¬ 
ously  exporting  the  entropy  it  produces  to  its  environ¬ 
ment.  As  a  result,  the  entropy  of  the  structure  plus  the 
environment  increases  at  a  much  faster  rate  than  when 
the  system  was  still  on  the  thermodynamic  branch  of 
the  solution. 

The  spiral  (or  bar)  patterns  of  galaxies  have  many 
characteristics  of  a  typical  “dissipative  structure” .  First, 
as  we  have  shown,  a  quasi-stationary  spiral  mode  is 
maintained  by  the  opposing  effect  of  the  spontaneous 
growth  tendency  and  local  dissipation,  with  a  contin¬ 
uous  flux  of  energy,  angular  momentum  and  entropy 
through  the  system  carried  by  the  trailing  spiral  wave 
itself.  Second,  it  can  be  shown  that  the  formation  of 
spiral  structure  accelerates  the  speed  of  entropy  evo¬ 
lution  of  a  spiral  galaxy,  compared  to  that  of  a  uni¬ 
form  disk,  by  several  orders  of  magnitude  (Zhang  1992). 
Thirdly,  since  a  spiral  mode  is  a  global  instability  in 
the  underlying  basic  state  of  the  disk,  the  spontaneous 
emergence  of  the  spiral  pattern  (which  is  obviously 
a  global  symmetry-breaking  process)  happens  as  long 
as  the  disk  satisfies  certain  far-from-equilibrium  con¬ 
straints  (i.e.  the  basic  state  characteristics  must  allow 
the  linear  growth  rate  of  a  spiral  mode  to  be  greater 
than  zero).  Lastly,  the  characteristics  of  the  quasi- 
stationary  spiral  pattern  formed  are  determined  solely 
by  the  properties  of  the  basic  state,  and  not  by  the  ac¬ 
cidentals  of  the  external  perturbations.  This  last  point 
is  reinforced  by  the  N-body  simulations  of  tidal  spiral 
patterns  in  slightly  unstable  disks,  where  it  was  found 
that  after  the  initial  transient  state,  the  characteris¬ 
tics  of  the  tidally-induced  patterns  correlate  strongly 
with  the  properties  of  the  basic  state,  rather  than  with 
the  nature  of  the  encounter  (Donner  &  Thomasson 
1994).  These  characteristics  of  the  spiral  structure 
clearly  identify  it  as  an  example  of  a  “dissipative  struc¬ 
ture”  defined  by  Glansdorff  &  Prigogine  (1971),  and 
the  spontaneous  formation  and  stabilization  of  a  large- 
scale  spiral  mode  as  an  example  of  a  nonequilibrium 
phase  transition. 

III.  COMPARISON  OF  THE  HIERARCHI¬ 
CAL  CLUSTERING  AND  THE  SECU¬ 
LAR  EVOLUTION  PREDICTIONS 

Currently  the  working  paradigm  for  the  formation 
and  evolution  of  galaxies  and  structures  is  the  hierar- 
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chical  clustering  or  cold  dark  matter  (CDM)  model.  It 
has  demonstrated  successes  in  reproducing  the  angular 
spectrum  of  the  cosmic  microwave  background  (CMB) 
radiation  as  well  as  many  of  the  aspects  of  the  distri¬ 
bution  of  large  scale  structure  (see,  e.g.,  Bahcall  et  al. 
1999  and  the  references  therein).  However,  at  the  in¬ 
dividual  galaxy  level,  the  CDM  model  has  encountered 
serious  challenges  in  attempting  to  reproduce  the  ob¬ 
served  galaxy  properties.  The  standard  CDM  is  now 
replaced  by  the  ACDM  paradigm,  though  many  of  the 
problems  still  remain. 

In  a  recent  article,  Peebles  (2002a)  compared  the 
current  state  of  cosmology  with  the  state  of  physics 
at  the  turn  of  the  19th/20th  century,  and  commented 
that  several  known  problems  of  the  CDM  could  po¬ 
tentially  turn  out  to  be  the  same  type  of  “Kelvin-level 
clouds”  which  a  century  ago  resulted  in  the  revolution 
of  modern  physics,  i.e.  the  creation  of  relativity  and 
quantum  mechanics  theories.  These  problems  include 
“the  prediction  that  elliptical  galaxies  form  by  mergers 
at  modest  redshifts,  which  seems  to  be  at  odds  with  the 
observation  of  massive  quasars  at  z  ~  6;  the  prediction 
of  appreciable  debris  in  the  voids  defined  by  L*  galax¬ 
ies,  which  seems  to  be  at  odds  with  the  observation 
that  dwarf,  irregular,  and  L„  galaxies  share  quite  sim¬ 
ilar  distributions;  and  the  prediction  of  cusp-like  dark 
matter  cores  in  low  surface  brightness  galaxies,  which 
is  at  odds  with  what  is  observed”  (Peebles  2002a). 

Historically,  CDM  type  of  theories  were  invented 
partly  to  get  around  the  problem  that  there  does 
not  seem  to  be  sufficient  time  for  the  seeds  of  the 
anisotropies  observed  on  the  cosmic  microwave  back¬ 
ground,  about  one  part  in  105,  to  grow  into  the  non¬ 
linear  structures  we  see  today  by  gravitational  means 
alone,  which  requires  seeds  of  one  part  in  ~  103  at  the 
time  of  recombination  (z=1000).  Furthermore,  the  Big 
Bang  nucleosynthesis  model  also  requires  a  significant 
amount  of  non-baryonic  dark  matter  (Primack  1999)  if 
the  universe  is  flat  as  the  inflation  scenario  suggests. 

Given  the  partially  ad  hoc  nature  of  the  introduction 
of  the  CDM  (especially  since  after  30  years  of  search,  no 
evidence  of  the  existence  of  the  CDM  material  has  been 
found),  it  should  not  come  as  a  total  surprise  that  prob¬ 
lems  surface  when  observational  data  become  available 
to  allow  a  detailed  comparison  with  the  predictions  of 
CDM  model.  In  fact,  most  of  the  problems  of  the  CDM 
scenario  (which  Peebles  had  quoted  three  above,  and 
which  we  will  list  several  more  in  the  following)  can 
be  characterized  by  that  it  prescribes  a  medium  for 
structure  formation  which  is  too  clumpy  (or  can  easily 
become  too  clumpy)  on  small  scales,  yet  too  smooth 
on  large  scales.  For  example,  the  cusp  problem  and 
the  satellite  abundance  problem  are  all  both  due  to  the 
over-clumpiness  of  the  medium  on  small  scales,  so  is 
the  problem  of  small  disks  or  rapid  angular  momentum 
loss  during  disk  formation  (White  &  Frenk  1991):  on 
the  other  hand,  the  over-smoothness  of  the  medium  on 
large  scales  underlies  the  problem  of  its  inability  to  ac¬ 
count  for  early  quasar  formation,  the  early  formation  of 


giant  high  redshift  clusters  (Francis  et  al.  1997:  Steidel 
et  al.  1998;  Williger  et  al.  2002),  as  well  as  the  problem 
of  accounting  for  the  observed  bubble  and  void  appear¬ 
ance  of  large-scale  structure  (Geller,  &  Huchra  1988). 
Furthermore,  in  a  purely  bottom-up  structure  forma¬ 
tion  scenario  such  as  the  CDM,  it  is  very  difficult  to 
account  for  the  alignment  of  the  spin  axis  of  the  bright 
galaxies  in  a  cluster  (Ozernoy  1974a  and  the  references 
therein;  Kim  2002;  Plionis  et  al  2003),  the  alignment  of 
substructures  in  a  cluster  (Plionis  &  Basilakos  2002),  as 
■well  as  the  observed  galaxy-cluster-supercluster  align¬ 
ment  effect  on  the  largest  scale  (West  2001). 

In  what  follows  we  contrast  a  number  of  the  major 
predictions  of  the  hierarchical  clustering  (CDM) /merger 
scenario  with  that  of  the  secular  evolution  (SE),  focus¬ 
ing  on  individual  galaxy  properties,  and  compare  both 
predictions  with  the  known  observational  facts  when 
available.  Through  this  cross  comparison,  we  wish  to 
demonstrate  that  secular  evolution  is  indeed  a  much 
more  natural  paradigm  in  explaining  the  properties  (as 
well  the  evolution  of  these  properties)  of  the  observed 
galaxies: 

•  The  CDM  model  predicts  that  the  total  number  of 
galaxies  of  all  Hubble  types  per  comoving  volume 
should  decrease  with  time  due  to  merger  events; 
whereas  SE  predicts  that  the  comoving  number 
density  of  all  Hubble  types  should  remain  nearly 
constant,  and  the  number  counts  for  individual 
Hubble  types  should  evolve  according  to  the  mor¬ 
phological  transformation  picture.  Recent  large 
surveys  such  as  the  Caltech  Faint  Galaxy  Red- 
shift  Survey  (CFGRS)  and  Slone  Digical  Sky  Sur¬ 
vey  (SDSS)  have  shown  that  there  is  essentially  no 
evolution  of  the  total  number  density  of  galaxies 
per  co-moving  volume  between  z=0  and  1  (Cohen 
2002;  Yasuda  et  al.  2001). 

•  According  to  the  CDM  paradigm  the  field  galax¬ 
ies  should  evolve  faster  through  merger  process 
(since  merger  is  known  to  happen  more  frequently 
in  the  fields),  and  rarely  happens  in  clusters  due  to 
the  relatively  high  speed  nature  of  cluster  galaxy 
encounters.  SE  on  the  other  hand  predicts  that 
cluster  galaxies  should  show  a  faster  secular  evo¬ 
lution  rate  than  field  galaxies  due  to  the  tidal- 
interaction-enhanced  and  spiral-mediated  mass  re¬ 
distribution.  The  well-known  Butcher  Oemler  ef¬ 
fect  for  duster  galaxies  (Butcher  &  Oemler  1978a, b) 
as  well  as  the  observation  of  field  galaxies  (El¬ 
lis  1997)  indicate  that  the  cluster  galaxies  have 
a  much  large  morphological  evolution  rate  than 
field  galaxies.  This  enhanced  secular  evolution 
rate  for  cluster  galaxies  is  also  consistent  with  the 
observation  that  cluster  core  spirals  are  found  to 
have  more  noticeable  falling  rotation  curves  (Ru¬ 
bin,  Ford,  &  Whitmore  1988). 

•  In  the  CDM  paradigm,  there  is  no  direct  relation 
between  the  kinematics  and  energetics  of  the  stars 
and  gas  in  a  given  spiral  disk;  whereas  the  SE  the- 
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ory  predicts  that  stars  and  gas  should  appear  on 
the  same  energetic  hierarchy  at  the  1  kpc  spa¬ 
tial  scale  (Zhang  2002),  because  the  interstellar 
medium  receives  similar  amount  of  energy  injec¬ 
tion  per  unit  mass  from  the  spiral  density  wave  as 
the  stellar  component.  The  observations  indicate 
a  clear  correlation  of  stellar  and  gaseous  kinemat¬ 
ics  on  1  kpc  spatial  scale  (Larson  1979,  1981;  Fleck 
1982). 

•  The  CDM  galaxy  formation  model  (Kauffmann 
1996)  prescribes  that  the  collapse  of  gaseous  ma¬ 
terial  within  the  clustering  dark  matter  halo  pro¬ 
duces  disks,  mergers  of  nearly  equal  mass  disks 
produce  ellipticals,  and  ellipticals  subsequently 
grow  disks  if  left  undisturbed.  One  consequence 
of  this  prescription  is  that  late-type  spirals,  which 
have  a  large  disk-to-bulge  ratio,  should  have  older 
bulges  than  do  early-type  spirals,  since  to  have  a 
larger  disk  the  galaxy  must  have  been  undisturbed 
and  be  able  to  accrete  gas  for  a  longer  time,  con¬ 
trary  to  the  observed  trend  (Wyse  et  al.  1997); 
the  SE  theory  on  the  other  hand  predicts  that 
the  early  type  bulges  are  older  since  they  have  ac¬ 
creted  mass  for  a  longer  time  and  most  of  the  ac¬ 
creted  mass  is  in  stellar  form  from  the  local  vicin¬ 
ity,  since  the  central  region  formed  stars  early  due 
to  the  deeper  potential  well  and  shorter  dynamical 
time. 

•  The  CDM  theory  does  not  predict  any  correla¬ 
tion  of  the  disk  heating  time  scale  and  the  angu¬ 
lar  momentum  transport  time  scale;  whereas  in 
the  SE  scenario  these  two  time  scales  are  found 
to  be  tightly  correlated  (see  equations  2  and  4). 
The  observed  continuous  change  of  bulge  to  disk 
ratio  and  other  type  of  bulge-disk  connections 
(Courteau  et  al.  1996)  attest  to  the  correlation 
of  heating  and  angular  momentum  transport  time 
scales  during  the  morphological  evolution  of  the 
disks. 

•  The  disk  galaxies  produced  in  the  CDM  simula¬ 
tions  are  much  too  small  and  rotate  much  too  fast 
compared  to  the  observed  galaxy  disks  (White  & 
Frenk  1991;  Navarro  &  Steinmetz  2000).  This 
problem  is  made  especially  acute  by  the  observa¬ 
tional  fact  that  the  observed  damped  La  systems 
at  high  z  (Wolfe  2001),  which  are  the  most  likely 
candidates  for  primordial  disks,  are  usually  quite 
large.  SE  starts  from  the  initial  condition  of  large 
disks  at  high  z,  which  become  the  progenitors  of 
today’s  large  early  type  galaxies. 

•  The  halo  merging  scenario  tends  to  create  a  cuspy 
mass  distribution  in  the  center  of  a  proto  galaxy. 
The  observed  young  galaxy  candidates  in  the 
nearby  universe,  i.e,.  the  low  surface  brightness 
(LSB)  galaxies  (Impey  &  Bothun  1997)  and  dwarf 
galaxies  are  found  not  to  possess  such  cuspy  cores 
(de  Blok,  McGaugh,  &  Rubin  2001).  Even  invok¬ 
ing  maximum  feedback  does  not  solve  the  cusp 


problem  (Gnedin  &  Zhao  2002).  The  SE  scenario 
on  the  other  hand  starts  from  the  more  flattened 
morphological  distribution  of  a  LSB  disk,  and  the 
central  density  of  a  galaxy  is  increased  gradually 
through  a  Hubble  time  of  evolution. 

•  The  CDM  simulations  of  galaxy  formation  have 
difficulty  fitting  simultaneously  the  observed  zero 
point  of  the  Tully-Fisher  (TF)  relation  and  the 
local  luminosity  function,  a  problem  related  to 
the  small-disk  problem  mentioned  above  (White  & 
Frenk  1991).  Furthermore,  the  CDM  model  natu¬ 
rally  predicts  that  the  luminosity  of  an  individual 
galaxy  is  proportional  to  the  third  power  of  its 
circular  velocity  (White  1997;  Dalcanton,  Spergel 
&  Summers  1997);  To  arrive  at  the  observed 
fourth  power  TF  relation  requires  fine-tuning  of 
the  feedback  and  cooling  parameters  (van  den 
Bosch  2000).  The  SE  scenario  on  the  other  hand 
predicts  that  the  entire  Hubble  sequence  from 
the  spiral  disks  to  disky  ellipticals  should  follow 
roughly  the  same  Tully-Fisher /Faber  Jackson  re¬ 
lation  (Zhang  1999).  No  modified  Newtonian  dy¬ 
namics  is  needed  to  explain  the  fact  that  LSB 
disks  fall  onto  the  same  TF  relation  as  normal 
spirals,  as  long  as  the  decrease  in  surface  bright¬ 
ness  of  a  LSB  galaxy  is  compensated  by  the  in¬ 
crease  of  dynamical  mass-to-light  ratio,  as  has 
been  found  to  be  the  case  for  low  surface  bright¬ 
ness  galaxies  (Zwaan  et  al.  1995).  Similar  trend 
also  holds  for  early  type  galaxies  the  increase  in 
surface  brightness  in  offset  by  the  decrease  in  dy¬ 
namical  mass-to-light  ratio  (Broeils,  Haynes,  & 
Baumgardt  1993). 

•  In  the  CDM  paradigm,  large  galaxies  are  formed 
out  of  the  mergers  of  smaller  ones,  and  therefore 
should  form  last.  Recent  studies  of  high  z  galaxies 
have  shown  that  exactly  the  opposite  is  observed: 
i.e.,  the  larger  the  mass  of  a  galaxy,  the  shorter  the 
time  scale  of  its  formation  (Thomas,  Maraston  & 
Bender  2002;  Boissier  et  al.  2001).  The  stellar 
population  in  the  lower  mass  disky  ellipticals  are 
found  to  be  systematically  younger  and  have  more 
extended  history  of  star  formation  than  that  in  the 
high  luminosity  boxy  ellipticals  (Bender  &  Saglia 
1999).  In  the  SE  scenario  large  galaxies  tend  to 
form  quicker  because  of  the  more  rapid  gravita¬ 
tional  collapse  to  form  disks  and  the  faster  rate  of 
evolution  due  to  a  more  prominent  spiral  structure 
on  the  massive  disk. 

•  The  observed  elliptical  galaxies  come  in  two  types. 
While  the  boxy  giant  ellipticals  have  characteris¬ 
tics  which  indicate  that  they  are  likely  to  be  the 
product  of  mergers,  i.e.,  being  pressure  supported, 
having  multiple  nuclei,  having  cuspy  cores  in  their 
central  density  profile,  having  two  populations  of 
globular  clusters  of  different  colors,  and  are  radio 
and  X-ray  loud,  etc.,  the  more  numerous  diskv 
type  ellipticals  on  the  other  hand  are  mostly  rota- 
tionally  supported,  having  power-law  central  den- 
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sity  profiles,  show  little  evidence  of  mergers,  hav¬ 
ing  only  one  globular  cluster  population,  and  are 
radio  and  X-day  quiet  (Bender  et  al.  1989;  Faber 
et  al.  1997;  Lee  2002).  It  is  our  belief  that  this 
dichotomy  of  characteristics  is  due  to  that  boxy 
ellipticals  are  the  true  merger  products,  whereas 
disky  ellipticals  are  produced  mostly  by  secular 
evolution, 

•  N-body  simulations  have  shown  that  stellar  merg¬ 
ers  tend  to  flatten  out  the  abundance  gradient 
from  that  of  the  observed  power  law  shape  (Mihos 
&  Hernquist  1994a).  Gas  rich  mergers  (i.e. ,  10% 
gas)  tend  to  create  a  distinctive  dense  core  (Mihos 
&  Hernquist  1994b)  which  is  once  again  not  ob¬ 
served  in  normal  ellipticals.  The  multi-component 
merger  scenario  (Weil  &  Hernquist  1994,  1996) 
is  not  relevant  to  the  explanation  of  the  gradual 
decrease  of  disk  galaxy  population  with  decreas¬ 
ing  redshift.  The  observed  density  and  abundance 
profiles  can  be  naturally  produced  and  enhanced 
through  secular  evolution  process. 

The  above  comparison,  as  well  as  the  fact  (see,  e.g. 
Peebles  2002b  and  the  references  therein)  that  it  ap¬ 
pears  difficult  to  reproduce  simultaneously  the  spec¬ 
trum  of  La  forest  (which  requires  plenty  of  small- 
scale  CDM  power)  and  the  individual  galaxy  properties 
(which  are  troubled  by  too  much  small  scale  power) 
indicate  that  the  true  underlying  structure  formation 
theory  may  not  be  purely  gravitational,  but  in  addition 
may  include  other  nonlinear  processes  such  as  primor¬ 
dial  turbulence  (von  Weisacker  1951;  Gamov  1952;  Oz- 
ernoy  1974b).  The  supersonic  shocks  associated  with 
turbulence  can  solve  the  time-of-growth  problem  of  the 
high-z  quasars  and  high  z  ellipticals  without  the  need 
of  dark  matter  (though  it  does  not  necessarily  exclude 
it);  it  also  naturally  produces  the  bubble  and  void  ap¬ 
pearance  of  the  large  scale  distribution.  Turbulence 
cascade  can  produces  the  observed  scale-invariant  mass 
spectrum  (Ozernoy  1974b)  just  as  the  CDM  theories. 
It  also  allow  small  scale  structure  to  be  present  but  not 
condense  and  collapse  to  form  cusps  and  cores.  The 
problem  with  the  primordial  turbulence  scenario  is  of 
course  why  such  supersonic  motion  has  not  left  signif¬ 
icant  imprints  on  the  cosmic  microwave  background: 
the  observed  CMB  fluctuation  of  10-5  does  not  leave 
room  for  significant  velocity  fluctuation  at  the  recom¬ 
bination  time  if  the  turbulent  matter  was  coupled  to 
radiation  then.  Recently,  the  first  evidence  of  primor¬ 
dial  turbulence  in  the  form  of  a  Kolmogorov  scaling 
relation  between  temperature  increment  and  angular 
separation  on  the  CMB  has  been  detected  (Bershad¬ 
skii  &  Screenivasan  2002).  The  issue  of  whether  this 
indicates  a  real  turbulent  state  of  the  matter  at  decou¬ 
pling,  or  else  it  is  rather  the  fossil  of  an  earlier  turbulent 
stage  before  decoupling  (Gibson  2000)  is  not  yet  clear. 


IV.  FUTURE  RESEARCH 

The  analytically  derived  evolution  rates  (2)  (4)  have 
been  quantitatively  confirmed  in  N-body  simulations 
(Zhang  1998),  in  the  sense  that  the  evolution  rates 
actually  measured  in  these  simulations  are  found  to 
be  nearly  identical  to  that  obtained  by  calculating 
the  evolution  rates  using  the  parameters  of  the  spi¬ 
ral  patterns  formed  in  these  simulations  to  plug  into 
the  above  analytical  evolution  rate  equations.  Us¬ 
ing  the  observed  spiral  arm  amplitude  of  30-50%  (e.g. 
Elmegreen,  Elmegreen  &  Montenegro  1992)  and  pitch 
angle  of  about  20%  (e.g.  Drimmel  2000)  in  these  rate 
equations  leads  to  the  conclusion  that  significant  mor¬ 
phological  transformation  within  a  fraction  of  the  Hub¬ 
ble  time  can  be  obtained  from  these  rate  equation. 

However,  the  spiral  patterns  formed  in  the  past  2D 
simulations  such  as  in  Zhang  (1996,  1998)  are  often 
found  to  have  the  amplitude  range  of  10-15%,  which 
leads  to  slower  secular  evolution  rate.  One  reason  of 
course  is  that  these  simulations  are  carried  out  for  per¬ 
fectly  isolated  galaxies,  which  in  reality  are  rare,  and 
which  when  exist  are  also  observed  to  have  slower  evo¬ 
lution  rates  than  their  cluster  or  group  counterpart. 
However,  we  suspect  that  there  might  exist  also  an 
artifact  in  these  past  2D  simulations  of  spiral  galaxy 
evolution,  especially  when  the  adopted  (bulge+halo) 
mass  to  disk  mass  ratio  is  high  (a  comparison  of  the 
2D  N-body  simulation  results  of  Zhang  1996  and  Zhang 
1999  shows  clearly  the  effect  of  a  large  rigid  spheroidal 
component  on  the  spiral  pattern  longevity  and  on  disk 
mass  accretion  rate,  as  well  as  the  plain  but  not  always 
obvious  truth  that  one  cannot  have  long-lasting  spiral 
structure  and  fast  morphological  transformation  at  the 
same  time:  kind  of  like  one  cannot  both  have  the  cake 
and  eat  it!).  Recent  studies  by  Athanassoula  (2002) 
found  that  making  the  halo  active  leads  to  enhanced 
bar  formation.  It  is  thus  expected  that  an  enhanced 
spiral  formation  and  enhanced  secular  evolution  rate 
will  also  be  obtained  in  a  full  3D  and  live  halo  simu¬ 
lation  of  spiral  disks.  Placing  a  galaxy  under  the  tidal 
influence  of  neighboring  galaxies  in  realistic  group  or 
cluster  environment  obviously  will  also  help  to  obtain 
a  larger  spiral  amplitude  and  thus  a  higher  evolution 
rate. 

Another  issue  which  can  be  explored  by  a  3D  simu¬ 
lation  is  the  secular  evolution  caused  by  a  skewed  mass 
distribution  such  as  found  in  the  high-z  field  galaxies  in 
the  Hubble  Deep  Fields.  These  skewed  3D  mass  distri¬ 
bution  should  produce  the  same  kind  of  phase-shift  and 
torque  relations  as  in  spiral  galaxies.  The  secular  evo¬ 
lution  produced  by  these  structures  may  also  played 
a  role  for  the  direct  formation  of  some  high-z  ellipti¬ 
cal  galaxies  (Peebles  2002c  and  the  references  therein) 
without  going  through  the  disk  formation  phase. 

The  secular  evolution  theory  we  presented  here  de¬ 
scribes  a  morphological  transformation  process  of  in¬ 
dividual  galaxies.  By  itself  it  does  not  uniquely  spec¬ 
ify  a  cosmology,  though  it  does  hint  at  the  elements 
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and  consequences  of  such  a  cosmology.  For  example, 
it  favors  large  disks  to  form  early,  and  subsequently 
undergo  morphological  transformation  mainly  due  to 
the  mediation  of  global  structures  such  as  spirals,  bars, 
and  three  dimensional  twisted  isophotes.  Environmen¬ 
tal  effect  accelerates  the  evolution  speed,  but  it  op¬ 
erates  mainly  through  the  mediation  of  internal  global 
structure,  and  not  through  the  actual  merging  of  galax¬ 
ies.  As  we  gradually  eliminate  uncertain  elements  from 
our  knowledge  of  the  galaxy  formation  process  we  will 
be  able  to  better  constrain  the  elements  of  earlier  cos¬ 
mological  processes. 

I  would  like  to  express  my  sincere  gratitude  to  pro¬ 
fessor  Hong  Bae  Ann,  professor  Hyung  Mok  Lee,  as  well 
as  the  Local  Organizing  Committee  for  inviting  me  to 
participate  in  this  fruitful  workshop  which  I  benefited 
much  from.  The  writing  of  the  current  manuscript  is 
supported  in  part  by  funding  from  the  Office  of  Naval 
Research. 
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